O-TT transition in magnetic triplet superconductor Josephson junctions 
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We examine a Joseplison junction involving two arbitrary equal-spin-pairing unitary triplet su- 
perconductors and a ferromagnetic tunneling barrier. Using perturbation theory, we show how the 
interaction of the barrier moment with the spin of the tunneling triplet Cooper pairs can reverse the 
sign of the Josephson charge current. This also results in a Josephson spin current, which contains 
a phase-independent contribution due to reflection processes at the barrier. We verify our analytic 
predictions using a non-perturbative Bogoliubov-de Gennes method. 
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Introduction. The interplay of superconductivity and 
magnetism is an enduring enigma of condensed matter 
physics. Over the last 25 years, many fascinating in- 
sights into this problem have been made in the study 
of singlet superconductor (SC) Josephson junctions with 
ferromagnetic (FM) tunneling barriers (H- For example, 
as the barrier width is increased, the usual Josephson 
current Ij vs phase relationship Ij — \ sin(0) becomes 
I J = I^ttI sin((/) + tt). This so-called O-tt transition is evi- 
dence of oscillations of the singlet SC correlations in the 
tunneling region 0, S]- A remarkable feature of such 
junctions is the presence of triplet SC (TSC) correlations 
induced by a proximity effect [J, 0j 01 > with the realized 
triplet pairing states dictated by the details of the FM 
barrier and the bulk SCs. Due to the likely intimate 
connection between triplet superconductivity and mag- 
netism , it is interesting to consider the case where the 
TSC pairing state can be chosen independently of the 
FM barrier. Despite the growing interest 0, [1, H 
in TSC Josephson junctions prompted by the discovery 
of Sr2Ru04 [lH, the study of such TSC-FM-TSC (TFT) 
junctions is still in its infancy. Recently, a novel O-tt tran- 
sition in a specific TFT junction was predicted [1, 
where the dependence of Ij upon the orientation of the 
FM moment indicates that it couples to the spin of the 
tunneling Cooper pairs. 

In this letter, we use perturbation theory [l^ to obtain 
the Josephson charge current through a TFT junction for 
arbitrary choice of unitary equal-spin-pairing TSCs. We 
predict that the O-tt transition found in Ref.s|3 and0 is 
always present for sufficiently large magnetization, and is 
due to the spin-flipping of tunneling triplet Cooper pairs. 
This also produces a Josephson spin current which 
has opposite sign on either side of the junction and a 
novel phase-independent contribution due to reflection 
processes. A non-perturbative Bogoliubov-de Gennes 
theory is used to demonstrate the universal character 
of our predictions, and that resonant tunneling through 
an Andreev bound state (ABS) does not qualitatively 
change the understanding of the O-tt transition [l^ [I3| ■ 

Perturbation theory. The Hamiltonian describing the 



TFT Josephson junction is written H = 'Hi,+T-LT^+Ht^,^ + 
Hra!- The terms and Hr respectively describe the bulk 
TSCs on the left and right side of the barrier: 



' ^..k (1) 



where V,.,k = {Cu,k,V^'^,KV^l,KV^i,KO'^ ^'^'^ ^U,-^ 
i^u k a) fermion creation (annihilation) operators, Cj^^k 
is the bare dispersion in the I'-hand TSC, and d^ k = 
Ajy kX are the triplet order parameters of the two TSCs. 
Both TSCs are in an equal spin-pairing state with respect 
to the z-axis and are unitary (i.e. the triplet conden- 
sate has no net spin) [l^. The gap in each spin sector 
is A^.,k.<T = -f7|A^,k|e*('^-+^-''^ where (j)^ is the global 
phase of the z^-hand TSC and 9^^]^ is an internal phase 
specifying the pairing state, obeying 6'i/^_k = Oi^.k + T^- As 
our results depend only on the spin state of the Cooper 
pairs, any variation of the orbital part of the gaps near 
the barrier will not qualitatively alter our conclusions. 

The two TSCs on each side of the barrier are linked by 
the tunneling Hamiltonian 



-i',k,(T''i',k' 
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where the subscript — i/ — R{L) when v = L{R). For a 
magnetically-active barrier, we must also include reflec- 
tion prcoesses fl7| : 
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Althogh reflection processes do not contribute to the 
charge current, spin- flip reflection may contribute to a 
Josephson spin current, as the spin- flip of a reflected 
Cooper pair changes the total spin in the TSC by ±2?i. 

In general, the matrix elements for spin-preserving tun- 
neling Tj^kk'' spin- flip tunneling T'j^k'k'; ^^'^ spin- flip 
reflection i?^'kk' ^''^ different. It is possible to derive 
expressions for the matrix elements from a more funda- 
mental Hamiltonian 17|, but here we will motivate a 
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phenomenological form. By Fermi's golden rule we have 
- and TZ-''- ^ \K:^;L'\' in the tun- 

nehng limit, where T"'"^' < 1 and W'^'^ < 1 are the 
transmissivity and spin-flip reflectivity of the barrier re- 
spectively. We also require that a tunneling or reflected 
quasiparticle acquires the same phase as in the exact so- 
lution. Following Ref. we consider the example of a 
purely FM barrier of (5-function width (appropriate for 
an atomically-thin barrier) at z = 0. We assume that 
the FM barrier moment M lies in the x-y plane at an 
angle r/ to the a;-axis. We hence use the Ansatz 



i/,k,k' 



(4) 
(5) 
(6) 



where M = g^.B\M.\/h^ 



with \ 



locity along the (OOl)-direction in the z/-hand TSC, j-sp, 



the Fermi ve- 

Ts, 

Tgf and Rsf are real constants and the 0{±kzk'^) guaran- 
tees that the transmitted or reflected quasiparticle moves 
away from the barrier The k-dependence of T^p, T^f 
and Rsf is irrelevant for our argument and is neglected. 
As rotating the spin coordinates about the x axis leaves 
the TSCs unchanged, our results for the charge current 
hold for any moment making an angle rj with the d^^k 
vectors. The spin current results also hold, but with 
corresponding rotation of the polarization. A schematic 
diagram of the junction is shown in Fig. ((T])(a). 

We define particle currents in the two spin sec- 
tors of each TSC by Ii,_a — ~i^{dtN„^a{t)) where 

N^A^) = EkcI,k,a(*)Cz.,k,a(0 and v = L{R) as a sub- 
script imphes v = — 1(-|-1) elsewhere. We calculate 
by expanding the S matrix to lowest order in Htun + '^,-ct, 
hence treating the tunneling and reflection processes as a 
perturbation of Ho = Wl + [I4I, which is justified for 
small T^'^-^i and -R^'kk'- "^^^ Kubo formula then gives 

- -i^y Jl^ dt'{[dtN,^c.{t),n,ut') + nut')])- 



interaction picture, we 

where we define the oper- 
t and 



k,.Wc,,k',.'(0 
i.,k,^WCi,,k', <;'(*)■ We consider 



Working within the 
write dtNy a{t) — 

ators Al^'-'it) = ^^.k 
B'/it) = Ek.k'^^;f " 

only the zero bias case, where the current through the 
spin-a sector of the i/-hand TSC is due only to the 
Josephson effect and may be conveniently expressed 
as li,^ = 2vIuy{^1^^{Q) + -^1^{Q)]. The retarded 
correlation functions and ^"^^(0;) give the 

tunneling and reflection contributions respectively. 
After using Wick's theorem to expand the two-particle 
correlators in the corresponding Matsubara functions 



make the analytic continuation iojn 



ATrAl'<-(T)Al^^' m and 
/o'' dre*""" (r,S-"'"(r)B-"^"(0)), we 



dre" 



the retarded functions. Substituting Eq.s into the 

expressions for the retarded correlation functions at 
= 0, we obtain the particle current 



I'' = 



E 

k,k' 



M4 E_,,_uE, 



k£/iy,k' 



X sin((/) -f j/[6l^,k' -6'-i.,k]) 



k,k' 



k£/iy,k' 



sin((/) -I- 2vari + v[9i,_\^i ~ 6'_iy^k]) 



k,k' 



Rlf |A^,kAt.,k' 
M2 ^,,k£;.,k' 



Fu,v,]^M'5^.:]^'9{-kzk'z) 



X sin(2m/ - [6l^,k - du^'V) 



(7) 



where Ejj\^ — a/ {^v,k — A*)^ + l^i^.kp is the excitation 
spectrum in the z/-hand TSC, (j) — (pB, ~ 4'l and 
FuyMM' = I]±["F{±-£^i^,k) - nF{E^\\^')]/[E^^^ =F £^i/',k'] 
with npiE) the Fermi distribution function. 

Eq. ^ is our first important result, as it contains all 
contributions to the current. The first term describes 
spin-preserving tunneling, where the Cooper pairs pre- 
serve their spin during the tunneling event, giving the 
usual Josephson result. The second term describes spin- 
flip tunneling, where the spin of the Cooper pair is re- 
versed by the coupling to the FM moment. Relative to 
spin-preserving tunneling, these Cooper pairs acquire a 
phase shift of 2ai'ri due to the spin-flip itself, and a fur- 
ther TT-shift arising from the intrinsic phase difference 
between the spin-t and spin-| condensates in the TSC. 
Lastly, we have the current due to Cooper pairs under- 
going a spin-flip when they are reflected at the tunneling 
barrier. As such, this term is independent of the TSC 
on the other side of the barrier, depending only upon the 
phase due to the spin-flip itself and the gap experienced 
by the reflected Cooper pairs. 

Charge current. From Eq. ([7]) we obtain the Josephson 
charge current Ij = —e{I^.^ + I;^ ^)- 



T o ( -^sp n'^s/ \ l^fl.kAi k'l CI 

// —2e — T — cos(277) — ^- > ■ ■ — Fb r k k' 



(^ky .kji 0{kzk'z) sm{4) + OR^y, - Ol.^) 



X dk 



(8) 



■ UJ - 



' iO"*" to obtain 



The first term in brackets corresponds to the spin- 
preserving contribution, while the second term is due to 
the spin-flip tunneling. The 77-dependence of the latter is 
due to the extra phase shifts for spin-flip tunneling: ig- 
noring orbital effects, a spin-a Cooper pair incident from 
the LHS undergoing a spin-flip during tunneling expe- 
riences an effective phase difference — (f> + — 2ar] 
between the two TSCs. The spin-flip current vs phase 
relationships in each spin channel are hence shifted by 
±4?7 with respect to one another. The interference be- 
tween the two spin channels results in the modulation of 
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FIG. 1: (color online) (a) Schematic diagram of the TFT 
junction, (b-d) The critical line r\a between the and tt states 
for various values of M in the (b) py, (c) pz, and (d) pz + ipy 
junctions. The tt state is realized for < Ir;], |7r — ?7| < 77c, 
while the state occurs for r]^ < \ri\ < n — rjc- The M S> 1 
behaviour of rjc at high T reflects the dominance of spin-flip 
processes; at low T, details of the bulk TSCs are important. 



the total spin-flip current by cos(2ry); this is analogous to 
the effect of the spin-dependent phase shifts for tunneUng 
between TSCs with misaligned d- vectors 011, S El- 
Spin- flip tunneling dominates Ij for large M: in this 
case, relevant for a half-metaUic barrier, we find Ij oc 
cos(2?7), and so the current changes sign at 77 = 7r/4 as 
the moment is rotated about the z-axis, i.e. there is a O-tt 
transition controlled by the orientation of the moment. 
As this originates solely from the spin structure of the 
triplet Cooper pairs it is an universal feature of unitary 
equal-spin-pairing TFT junctions, our second important 
result. To test this prediction, we consider TFT junctions 
where both TSCs are made of the same material, for the 
three choices of Py, Pz, and Pz + ipy orbital pairing sym- 
metry (the latter of most relevance to Sr2Ru04 [11]). In 
a model with spatially constant TSC gaps, we can include 
tunnehng and reflection processes to all orders [isl. [itI [l8| 
by solving the Bogoliubov-de Gennes equations to ob- 
tain the ABS energies E±^i^{(j),r], M,T) at temperature 
T. The free energy of the junction is then given by 
F = -/3EkEp=±(|fc-l/fcF)log(2cosh(/3£;p,k/2)); we as- 
sume a 2D circular Fermi surface in the y-z plane. The 
ground state of the junction is found by numerically min- 
imizing F with respect to </). For each junction, the min- 
imum of F lies at 15!) = or tt; we find a O-tt transition 
when the global minimum shifts from one of these values 
to the other. In Fig. (II|)(b-d) we plot the critical an- 
gle r]c at which this occurs in each junction for fixed M; 
according to Eq. ([8]), in the tunneling limit ijc = tt/A. 

For T sufficiently close to the transition temperature 
Tc, rjc always approaches the tunneling limit results as 



M is increased. For the Py junction [Fig. (IT|)(b)], r]c 
shows only weak T-dependence at fixed M, consistent 
with Eq. ([5]). Inthepz and -I- junctions [Fig. ([T|)(c) 
and (d) respectively], however, 77c varies significantly with 
T. The key difference between the Py junction and the 
Pz and Pz + ipy junctions is the absence of a zero energy 
ABS in the former. In the latter, there is a zero energy 
ABS at = TT for any choice of tunneling barrier, which 
raises the free energy of the tt state, thereby suppressing 
the O-TT transition. The strongest deviations from Eq. ([5]) 
therefore occur for the Pz junction, as here a zero energy 
ABS forms for all k; for the Pz+iPy junction, in contrast, 
a zero energy ABS forms only when k || z. Our pertur- 
bation theory results are recovered at higher T in these 
junctions due to the suppression of multiple-Cooper pair 
tunneling processes, which are a key feature of tunneling 
through a zero energy ABS [H, [l9| . 

Spin current. The spin of a triplet Cooper pair also 
allows a zero-bias Josephson spin current to flow across 
the junction M, [^, ■ The spin current is polarized along 
the z-axis [10], and so from Eq. ^ we obtain the spin 
current in the j/-hand TSC I^-j = f (/^{^ - I;,^^): 



j-S,z _ 



xSk„,k'^9ikzk'^) cos{(t) + Oj. 



fi?,k — ^'L,k')-^fl.,L,k,k' 

+hL' sm 2?7 2^ — ^ — — 

x^k|i,kj|6'(-fc^fci)cos(6'^,k - (9) 

The first term is from spin-fiip tunneling, while the 
second (/)-independent term is due to spin-fiip refiec- 
tion [2^. The spin-dependent phase shifts of the spin- 
flipping Cooper pairs are responsible for driving the spin 
current, again in analogy to the spin current between 
TSCs with misaligned d vectors 0,13, 113] • This implies 
I^'j — for rj = n7r/2, n G Z, as the relative phase be- 
tween the spin-flip currents in each spin channel is then 
2mT. As shown in Fig. ([2]) (a), the spin current reverses 
sign across the barrier: the spin current carried by the 
tunneling Cooper pairs reverses on the spin-fiip, while the 
spin-flip reflected Cooper pairs in each TSC carry oppo- 
site spin current as they move in opposite directions. 

Eq. ([9|) may be simplified for the three junctions intro- 
duced above. By energy conservation we have k' = k 
in the tunneling term and k' = k = {kx,ky,—kz) in 
the refiection term. Furthermore, we set Tsf = Rsf as 
here T'^^~'^ — TZ"''~"'; the amplitude of the cosine term in 
each contribution to Eq. ^ is then identical. We hence 
find Ij^'j oc 7 -|- cos(0) where 7 is an orbital-dependent 
constant due to the phase shift AO^^y^ — 9i,_k — 9^"^ 
perienced by specularly-reflected Cooper pairs. For the 
Py junction we have Q^^ = ^j/,k, and so there is no ex- 
tra phase shift upon refiection, giving 7 = 1. For the 
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As J = for 77 = and r\ = tt/2, the different sign of Ij 
at these angles is however a robust equiUbrium feature. 

Conclusions. We have analyzed the Josephson cur- 
rents through a TFT junction for any choice of unitary 
equal-spin-pairing TSCs. We predict that the sign of 
the charge current is controlled by the relative impor- 
tance of spin-flip to spin-preserving tunneling. Spin-flip 
processes also produce a Josephson spin current, with a 
phase-independent term due to reflection. Our results 
reveal the importance of the Cooper pair spin as a novel 
degree of freedom in TSC Josephson junctions. 

We thank Y. Asano, J. Linder, D. K. Morr and M. 
Sigrist for useful discussions, with special thanks to B. 
Rosenow and C. Timm. 



FIG. 2: (color online) (a) Cartoon of spin-flip tunneling (solid 
line) and reflection (broken line) processes contributing to the 
spin current, (b-d) The 2-component of spin current on the 
LHS for various values of 7; at M = 2 and T — OATc in the 
(b) py, (c) pz, and (d) + ipy junctions. The legend in (b) 
is for all plots. Ao is the maximum gap magnitude at T = 0. 



Pz junction, in contrast, all reflected Cooper pairs ex- 
perience a TT phase shift and therefore 7 = — 1. In the 
Pz + iPy junction, A9^ \^ — n ~ 2 arctan(fcy/fcz) depends 
upon k; integrating across the Fermi surface we find 
— 1 < 7 < 0. We again verify these predictions within the 
Bogoliubov-de Gennes theory for spatially-constant TSC 
gaps. Solving for the scattering wavefunctions SHE!, 
we obtain the Andreev reflection amplitudes af^^^'^J^) for a 
spin-tj electron-like (hole-like) quasiparticle incident from 
the j^-hand-side Andreev-reflected as a spin-cr' hole-like 
(electron-like) quasiparticle. Following Ref. |2l|, we write 



the spin current in terms of the ^ 

\k 



h{he) 



tS,Z 



1 



^a{a^^,,,(k,zc.„) - a%^„{\^,iu:^)) (10) 



We verify the relation /^'j = —I^'^j (not shown), and also 
find excellent agreement with the tunneling Hamiltonian 
predictions for 7 in all three junctions, see Fig. ©(b- 
d). The role played by reflection processes in the spin 
transport is our third important result. 

We have not accounted for the transfer of spin to the 
barrier moment when the spin current is non-zero. This 
can be physically justified if the barrier is in contact with 
a spin reservoir, allowing the diffusion of the transferred 
spin. In the absence of such a reservoir, we speculate 
that the barrier moment will precess about the x-axis, as 
the spin current only has a spin polarization || d x M. 
This very interesting matter requires a non-equilibrium 
treatment, which is beyond the scope of the present work. 
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